Influence of the Barrier Shape on Resonant Activation 
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The escape of a Brownian particle over a dichotomously fluctuating barrier is investigated for 
various shapes of the barrier. The problem of resonant activation is revisited with the attention on 
the effect of the barrier shape on optimal value of the mean escape time in the system. The char- 
acteristic features of resonant behavior are analyzed for situations when the barrier switches either 
between different heights representing erection of a barrier and formation of a well, respectively, or 
it proceeds through "on" and "off" positions. 
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I. INTRODUCTION 



Nonequilibrium systems driven by noises are known 
to display plentitude of noise-induced phenomena ^. 
Among those, stochastic resonance is manifested when 
the response of a nonhnear system to a signal is en- 
hanced by the presence of noise. Another resonant phe- 
nomenon can be observed for thermally activated sur- 
mounting of a potential barrier with a randomly fluc- 
tuating shape. For certain values of characteristic pa- 
rameters of the noise, the transport over the barrier is 
facilitated, i.e., the mean escape time of a particle ex- 
hibits a minimum as a function of the parameters of the 
barrier fluctuations. Such a problem of thermally acti- 
vated escape within randomly fluctuating potentials oc- 
curs in a wide variety of contexts with examples ranging 
from molecular dissociation in strongly coupled chemical 
systems [Q or the model dynamics of the dye laser ^ to 
selective ion pumps in biological membranes Q|. In ev- 
ery case activation happens due to thermal fluctuations 
and after classical Kramers theory can be expressed 
by an Arrhenius dependence for the mean lifetime of the 
metastable state, W ^ exp{AG/kBT), where AG stands 
for the activation energy. The strong effect of the sur- 
roundings ( "environmental noises" ) can be readily under- 
stood since even small variations 5G in AG will greatly 
affect W provided SG > ksT. For a complex nonequi- 
librium system the potential experienced by the Brow- 
nian particle cannot, in general, be regarded as static 
but rather, as influenced by random fluctuations whose 
characteristic time scale may be comparable with the du- 
ration of the diffusion over the barrier. Random effects of 
the environment can thus be envisioned as barrier alter- 
nating processes that can modulate the escape kinetics 
in the system. If the barrier fluctuates extremely slowly, 
the mean first passage time (MFPT) to the top of the 
barrier is dominated by those realizations for which the 



barrier starts in a higher position and thus becomes very 
long. The barrier is then essentially quasistatic through- 
out the process. At the other extreme, in the case of 
rapidly fluctuating barrier, the mean first passage time 
is determined by the "average barrier" . For some par- 
ticular correlation time of the barrier fluctuations, it can 
happen however, that the mean Jkinetic rate of the pro- 
cess exhibits an extremum |, 



10 that is a signature of 
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resonant tuning of the system in response to the external 
noise. 

In this communication we present numerical results for 
the mean first passage time over a fluctuating barrier for 
the model system with an arctan potential barrier subject 
to dichotomous fluctuations. Our model belongs to the 
class of "on-off" models discussed in a seminal paper by 
Doering and Gadoua and further analyzed by Bogufia 
et al. [jll] . A distinctive characteristics of this model is 
that part of the time the barrier is either switched off {i.e. 
it becomes flat), or the switching is performed between 
the barrier and a well, so that the particle can essentially 
roll rather than climb during these times. 

Our major goal was to examine variability of the mean 
first passage time in parameter space, that is, to deter- 
mine the position of its minimum (c/. Fig. |l|) and to 
study the dependence of depth and width of this mini- 
mum on the parameters describing the barrier shape and 
correlation time of the environmental noise. 

In Section II, a brief statement of the "archetypal" 
resonance activation problem is presented. Section III 
discusses solutions to the problem posed for the arctan 
potential, pointing out special symmetry of this function 
that determines equality of the mean first passage time 
for convex and concave slopes approximating between the 
triangular, piecewise-linear and piecewise-constant forms 
of the potential. Numerical results are obtained and fur- 
ther analyzed to asses how sharp and persistent the res- 
onant region of the mean first passage time can be as a 
function of the noise correlation time and steepness of the 
potential slope. We conclude with some final remarks in 
Section IV. 
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In the above equations time has dimension of 
[length]^ / energy due to a friction constant that has been 
"absorbed" in a time variable. We are assuming a re- 
flecting boundary at a; = and an absorbing boundary 
condition at a: = L 

P(L,a±,i) = (3) 




log,J(L'/T)T] 



(4) 



The initial condition 



P{x,a+,0)^P[x,a^,Q) = -5{x) (5) 

expresses equal choice to start with any of the two con- 
figurations of the barrier. The quantity of interest is the 
mean first passage time 



FIG. 1: The resonant t(7) line and its asymptotes at low and 
high external noise frequencies. 



II. GENERIC MODEL SYSTEM 

We have considered an overdamped Brownian parti- 
cle [0 moving in a potential field between absorbing and 
reflecting boundaries in the presence of noise that modu- 
lates the barrier height. The evolution of a state variable 
x{t) is described in terms of the Langevin equation 



dx 



— = -v'[x) + v^m 



-T4(a;) + V2re(t). 



(1) 



Here £^{t) is a Gaussian process with zero mean and 
correlation < £,{t)£^{s) >= 6{t — s) {i.e. the Gaussian 
white noise arising from the heat bath of temperature 
T), ri{t) stands for a Markovian dichotomous (not nec- 
essarily symmetric) noise taking on one of two possible 
values a± and prime means differentiation over x. The 
correlation function of the dichotomous process satisfies 
< r]{t)r]{t') >= e~^'^'*~*'l(a+ - a_)^/4 where 7 stands 
for the flipping frequency of the barrier fluctuations {i.e. 
^ represents the characteristic correlation time of the 
dichotomous noise). Both noises are assumed to be sta- 
tistically independent, i.e. < ^{t)ri{s) >— 0. Equivalent 
to eq. (|l|) is a set of the Fokker-Planck equations de- 
scribing evolution of the probability density of finding 
the particle at time t at a position x, subject to the force 
-V±{x) = ~V'{x) + a±g{x) 



dt J [P{x,a+,t) + P{x,a^,t)]dx 


= r+(0)+T_(0) (6) 

with T+ and r_ being MFPT for (+) and (— ) configu- 
rations, respectively. MFPTs r+ and r_ fulfill the set of 
backward Kolmogorov equations M, |fl 



~JT±{x) + JT^{x) 

dV±{x) dT±{x) _^^dV±(x) 



dx dx ' dx"^ 
with the boundary conditions (c/. eq. (^) and (|^)) 

T±{x)\x=0^^ T±{X%=L=^ 



(7) 



(8) 



Although the solution of (0) is usually unique Q, a 
closed, "ready to use" analytical formula for r can be 
obtained only for the simplest cases of the potentials 
(piecewise linear and piecewise constant). More complex 
cases, like even piecewise parabolic potential V± result in 
an intricate form of the solution to eq. (^. Other situ- 
ations require either use of approximation schemes 
perturbative approach |^ or direct numerical evaluation 
methods [|4| . In order to examine MFPT for various po- 
tentials, a modified program | [T5[ | applying general shoot- 
ing methods has been used with part of the mathematical 
software obtained from the Netlih library. 



III. SOLUTION AND RESULTS 



dtP{x,a±,t) = d^\V±{x)+Td.,\P{x,a±,t) 

- jP{x,a±,t)+-fP{x,a^,t) (2) 



For 7 = and V+{x) = V-{x) = V{x), eq. (|) de- 
scribes an overdamped Brownian particle moving in an 
external static potential (c/. eq. (|l|) with ri{t) =0). In 
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this situation a formula for the MFPT reads 



cxp 



T 



L L 



exp 



T 
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and does not depend on a constant part of a potential. 
Another non-trivial property in this case is equality of 
times 

r(C/,0 ^ L) = t(F,0 ^ i), (10) 

true for any two potentials V{x), U{x) that fulfill 

U{x) ^1-L-V{L~ x) (11) 

and are continuous functions mapping [0, L] — > [0, H]. 
The above property can be simply proven after noticing 
that the difference 



'(y,0 



L L 

L) - t{U, 0^L) = J I f{7i Odvd^ (12) 





with 



cxp 
T 



viv)-vm 



exp 



T 



(13) 



vanishes for any points (p, q) and (p' , q') such that (p', q') 
is the image of {p, q) in the reflection at the line rj ~ L — ^ 



f{p,q) = -f{p',q')- 



(14) 



To capture the features of the mean first passage time 
as a function of the shape of the barrier and character- 
istic correlation time of the barrier fluctuations, we have 
studied the Brownian diffusion problem in the arctan po- 
tential fluctuating randomly between two conflgurations 
V+{x) and V-ix) 



V±ix) 



i/±arctan(aa;/L + /3) 
arctan(Q! + (3) ~ arctan(/3) 

i?±arctan(/3) 
arctan(Q! + (3) — arctan(/3) 



(15) 



Since switching between both configurations is modelled 
by a Markovian two-state process, the potential expe- 
rienced by a particle remains in a given configuration 
for an exponentially distributed time, before flipping to 
the other state. The change in steepness and shift of 
the potential slope is controlled by parameters a and 
f3. By using this particular form of the potential, we 
are able to modulate the steepness of the slope when 




FIG. 2: Potentials of arctan type for: A f3 — —a, B /3 
-a/2, C P = 0. 



holding the same, constant values of V±{x) at the bot- 
tom of the potential (x — 0, V±(0) — 0) and at the top 
{x = L, V±{L) = H±) of the barrier (see Fig. |). More- 
over, as it is clearly seen from eq. (|l^), for /3 = and 
/3 — —a, the potential satisfies the property eq. ( [Ill ) with 



V±{x)\p^o^H±-V±{L~x)\p^. 



(16) 



Thus, under the same type of barrier flipping process, the 
MFPT results are identical for the convex (/3 = 0) and 
concave (/3 = —a) potential slopes. 

The standard numerical analysis (l|, |l§ of the MFPT 
has been performed for various forms of the potential 
eq. ( [l5| ) with the slope changing accordingly to a = 2% 
i = 0, . . . , 10 and J? = -a (Fig. | A), [3 = -a/2 (Fig. |B) 
and [3 = (Fig. || C). The height of the potential barrier 
H± has been switching between either two symmetric val- 
ues, H± = ±8r, or "on" and "off" barrier positions, i.e. 

= 8T, H- = 0. Numerical analysis reestablishes dis- 
cussed above equality of the MFPTs for potentials with 
(3 = —a and /3 = 0. It means that for this set of param- 
eters, the property eq. ( p^ observed for static potential, 
is also recovered for dichotomously fluctuating barriers 



r(C/,7)-T(V^,7) 



(17) 



Therefore, graphical presentation of the results relates 
arbitrarily to the convex potentials only {(3 = 0), having 
in mind that the same relationship t{^) is observed for 
the concave potentials with (3 = —a {cf. Fig. ||). 

Generally with increasing slope of the barrier, resonant 
frequencies are shifted to higher values suggesting that 
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more energy has to be pumped into the system in order 
to facihtate the escape of a Brownian particle over the 
fluctuating barrier. The minimal (resonant) values of r 
increase with a for V± — ±y potential and decrease for 
potentials switching between the "on-off" (V+ — V, V- — 
0) positions, thus documenting a different character of 
resonant activation in these two cases, as discussed by 
Bogufia et al. []lT| . 
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FIG. 4: The r{j) as a function of barrier fluctuation rate for 
arctan potential with H+ = 8T, H- = -8T, p = {p = -a). 



parameter a. 



FIG. 3: The r(7) as a function of barrier fluctuation rate for 
arctan potentials with H+ = 8T, H- = -8T, (3 = -a/2. 

For V± — ztV barriers asymptotic values of MFPT 
at 7 cxD do not change with the steepness: they re- 
main the same (Figs. ||, ^) and equal to as predicted 
by analytical results for piecewise linear and piecewise 
constant potentials j3[ ^ 0. However, they alter for 
the "on-off" potentials (c/. Figs. ^ and ^ displaying in- 
crease for f3 — —a/2 and decrease for /3 = (same for 
concave potentials /? = —a). The resonant line (Fig. of 
the t(7) dependence displays different features with zero 
and non-zero shift p. In particular, for /3 = the relative 
depths of minima measured from either the asymptotic r 
value at 7 = (ho), or from the r asymptote at 7 — > cxi 
(/loo) decrease with increasing a {cf. Figs. I and |). For 
non-zero (3 = —a/2 and both potentials, V± = ±V or 
V+ = V, V- = 0, the relative depth ho increases with a, 
as documented in Fig. ^ where the distance between the 
lines logio[r(7o)T/L2] and logio[T(7mm)T/-L2] is shown 
to increase with this parameter. At the same conditions, 
the depth of the resonant line /loo measured from the 
asymptote at 7 ^ 00 decreases for V± — ±V potentials 
and increases for the "on-off" potentials V+ — V, V- — 0. 
These findings are summarized in Figs. |^ and p| where log- 
arithmic values of 7™™, ri'jmin), 'r(7o), r(7oo) for the aU 
arctan type potentials are plotted versus varying slope 



H^=8T H=0 a=1 M-5 

H^=8T H =0 a=8 |l=-4 

H^=8T H =0 a=32 P=-16 

H=8T H =0 a=1024 |i=-512 




log,.[(L'/T)7l 



FIG. 5: The r(7) as a function of barrier fluctuation rate for 
arctan potentials with H+ = 8T, H- = 0, /3 = —a/2. 

Moreover, as can be easily seen in Figs, ^jand^, the res- 
onant width r/2 measured at hoo/2 {cf. Fig. |l|) increases 
significantly for the "on-off" potentials with increasing 
value of the slope parameter a. It is thus suggestive of 
an easier fine-tuning of the system to resonant conditions 
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O (IJ=-ct), H^=±8T 
l09,„[T„,„l-''/T] 
log,„[T(Y„,„)T/L=] 

log,„[T(Y„)T/L=] 

loa,„[^(T,.)T/L=] 




l09„[|LT)Tl 



FIG. 6: The r{-f) as a function of barrier fluctuation rate for 
arctan potential with H+ = 8T, H- = 0, P = {f3 = -a). 



FIG. 8: Location of minima, minimal and asymptotic values 
of the MFPT for arctan type potentials with /3 = (/3 = —a). 
Lines have been drawn to guide the eye. 



in this case by switching frequency of the barrier. The 
effect is also reproducible for /? = 0, although in this case 
the relative depth of MFPT minimum decreases, causing 
observable "flattening" of the resonant line. 



:-ct/2, H^=±8T 

- — iog,„[T„i„L°/r| 
« — log,„[T(T„„)T/L=] 

log,„[T(yJT/L'] 

log,„[T(Y_)T/L=] 



P=-0t/2, H =8T, H =0 

— ■ — log,„[Y„„L=/T] 

• log,„[T(Y„,„)T/L''] 

^ log,„[^(Y„)T/L=] 
' log,jT(Y )T/L'l 



FIG. 7: Location of minima, minimal and asymptotic values 
of the MFPT for arctan type potentials with /3 — —a/2. Lines 
have been drawn to guide the eye. 

The functional dependence t(7) has a typical charac- 
ter in all cases discussed in this work. In agreement with 



theoretical considerations ||, |T|, |T|, |T|, |T|, H, |T|' the 

asymptotic values t(7o) and t(7o^ for large and small 
correlation times of the barrier noise are recovered with 
a minimum located in between those two. For all slopes 
of barriers this characteristic behavior of the MFPT re- 
mains similar although the weakening of the effect, i.e. 
shift of the resonant "fmin to higher values and shallowing 
of the resonant line are observed with increasing steep- 
ness of the barrier. The effect of attenuation is stronger 
for V± — ±.V potentials with both, f3 — —a/2 or /3 = 
(/3 = —a). The resonant line flattens and slowly disap- 
pears with increasing a {i.e with increasing steepness of 
the barrier slope). The result is documented in Figs. |^ 
(upper panel) and ^, where the noticeable convergence 
between the values r(7min) and t(7oo) is observed for 
a = 1024. Numerical estimates of MFPTs for the al- 
ternating potential (V± ~ ±8T) yield T{"/min) ~ 0.3 and 
t(7oo) « 0.5 for a = 8, /3 -4 leading to h^. « 0.2. This 
value drops by an order of magnitude to hoc ~ 0.01 for 
a = 1024, /3 = —512, that corresponds to by two-order 
of magnitude higher steepness of the slope estimated as 
the derivative of the potential at point x = L/2. Thus, 
the result suggests a fairly robust character of the reso- 
nant activation that disappears, but only for the limiting 
step- function potentials. An apparent persistence of the 
phenomenon for V+ = 8T, V- = and P = -a/2 (Figs. | 
and 0, lower panel) is due to a different scenario of the 
barrier passage. In this particular case, an increasing a 
produces a wall at x = L/2 which has to be passed by 
a Brownian particle before being eventually absorbed at 
X = L. The particle experiences then a non-zero deter- 
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ministic force at this point, only. 

Potentials of the arctan type (eq. ( |l5| ) and Fig. ||) 
produce barriers approximating between triangular and 
rectangular shapes mostly used for analytical studies of 
the resonant activation effect |ll|, |l^, |2^. Although 
the resonant phenomenon has been detected and analyt- 
ically proven as a generic property in the case of dichoto- 
mously switching piecewise linear potentials |l^, it 
has been shown jT7| not to occur in the case of piecewise 
constant potentials, when the Brownian particle escapes 
over either the small or the large barrier at all times. 
In consequence of the latter, barrier crossing events and 
barrier fluctuations are two independent stochastic pro- 
cesses 0, m with the mean exit time decreasing 
monotonically from t Ki t-\. for 7 < to r « t_ 

for 7 > tZ^ ■ The qualitative behavior of r is then 
identical with the behavior of the exit time in so called 
kinetic models jl^, |9[ |2^; ^ which explicitly in- 
clude the escapes over both high and low barriers at a 
rate W± = cxp{—H±/T) and transitions between the 
two configurations of the potential. As it has been dis- 
cussed by Bier and Astumian , the kinetic description 
assumes an instantaneous adiabatic adjustment and ac- 
cordingly leads to the set of equations for the populations 
7T± of particles that feel the potential V± and have not 
yet escaped over the barrier 



dt 



TT±{t) ~ — W±'!T± — JTT± + 77rzp 



(18) 



The decay of the survival probability Tr{t) = 7r+ (t) +7r_ (t) 
is then determined by an average MFPT 



T - 



dt 



27+(W+ + VK-)/2 
W+W- + {W+ + W^_)7 

(19) 



which for slow barrier fluctuations (7 = 0) approaches 
the average of the MFPT for the alternative configura- 
tions and in the fast fluctuation limit (7 ^ 00 but re- 
specting 7 << T^^, where T7^ stands for the frequency 
of the escape attempts [|l3[ ]l^, |2^) coincides with the 
reciprocal of the average rate over the fluctuating bar- 
rier. For fluctuating rectangular potentials with a barrier 
placed at x — L/2, the closed expressions for the mean 
escape time Q for 7 = and 7 ^ 00 read 



-(70) 



T 
Z2 



(20) 



T 1 
r(7oo)-^ = 2 



T+T_ — 1 



(21) 



and are perfectly reproduced in numerical MFPT studies 
in the V± = ±y potential with parameters a = 1024, (3 = 
—a/2 approximating a step-like function 9{x — i/2). 
However, in the case of the "on-off" potential V+ = 
8T, = and the same set of parameters a, (3 numerical 
T value for 7 — > 00 differs from the eq. (^ij) by four order 
of magnitude (c/. Figs. ^ ^) displaying a persistent reso- 
nant activation vanishing slowly at much steeper barriers 
which eventually reach the limit of an infinite tangent at 
X = L/2. Similarly to the former results [ p^ , in the 
limiting case of a piecewise constant potential, the mean 
exit time approaches r « r_|_ for 7 < and r « r_ for 
7 > rJT^ displaying an inverse proportionality to the rate 
of the barrier fluctuations for intermediate rates. 



IV. SUMMARY 

In the foregoing sections we have considered the ther- 
mally activated process that occurs in a system coupled 
to an external noise source. The external stochastic pro- 
cess is responsible for fluctuations of the potential bar- 
rier which has been modelled by an arctan function of 
varying slope. As expected based on previous theoretical 
studies ^ , the phenomenon of resonant activation oc- 
curs typically in the system under broad circumstances 
of varying shape of the potential barrier and qualities 
of barrier fluctuations. In general, with increasing bar- 
rier steepness the resonance phenomenon becomes less 
"sharp" with the flat region of flipping frequencies around 
the resonant one. The time scale of the resonant activa- 
tion process is fairly insensitive to the potential param- 
eters except in the case of the "on-off" potentials, when 
broadening of the MFPT resonant line occurs suggest- 
ing most efficient tuning of the system to the resonance 
conditions. 
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